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1. Introduction 

Recently, we have obtained a suitable extension of the canonical fluctuation-dissipation 
relation involving the heat capacity C and energy fluctuations [1]: 

kBC = (3^5U^) + C{5l3JU), (1) 

which considers a system-surroundings equilibrium situation in which the inverse 
temperature (3^ = l/T^^ of a given thermostat exhibits non- vanishing correlated 
fluctuations with the total energy U of the system under study as a consequence of 
the underlying thermodynamic interaction, {5(3^5U) ^ 0. Clearly, Eq.([l]) differs from 
the canonical equilibrium situation due to the realistic possibility that the internal 
thermodynamical state of the thermostat can be affected by the presence of the system 
under study. This allows us to describe the fluctuating behavior of the system under 
more general equilibrium situations, rather than the ones associated with the known 
canonical and microcanonical ensembles. 

The fluctuation relation ([1]) possesses interesting connections with some challenging 
problems related to statistical mechanics, such as: (1) a compatibility with the existence 
of macrostates exhibiting negative heat capacities C < [H [2], a thermodynamic 
anomaly that appears in many physical contexts (ranging from small nuclear, atomic 
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and molecular clusters [21 SI El E] to the astrophysical systems [3 El El [ID]) associated 
with the existence of nonextensive properties [HI [121 [13 HI] ] (2) a direct application 
for the extension of available Monte Carlo methods based on the consideration of the 
Gibbs canonical ensemble in order to capture the presence of a regime with C < 
and avoid the incidence of the so-called super- critical slowing down [H [15] (a dynamical 
anomaly associated with the occurrence of discontinuous (first-order) phase transitions 
[16], which significantly reduces the efficiency of Monte Carlo methods [H]); (3) finally, 
a direct relationship with an uncertainty relation supporting the existence of some 
complementary character between thermodynamic quantities of energy and temperature 
[H HE] , an idea previously postulated by Bohr and Heisenberg [HI [20] with a long history 
in the literature [2T[ l22l [231 [2l] . 

Our aim in this work is to present a more complete study of the existing connections 
of the fiuctuation-dissipation relation ([T]). The core of our analysis is focussed on certain 
geometric aspects relating the present approach with other geometric formulations of 
fluctuation theory [25]. Such ideas straightforwardly lead to a geometric generalization 
of the Gibbs canonical ensemble describing a special family of equilibrium distributions 
recently proposed in the literature [261 [2Z], which can also be obtained from some 
known formulations of statistical mechanics, such as Jaynes' reinterpretation in terms 
of information theory [28], as well as Mandelbrot's approach based on inference theory 
[29]. 

2. A brief review 

2.1. Compatibility with negative heat capacities 

Our main motivation in deriving the fluctuation-dissipation relation was to arrive 
at a suitable extension of the known fluctuation relation: 



that is compatible with the existence of macrostates with negative heat capacities [Tl[2]. 
As discussed in many standard textbooks on statistical mechanics [16], the latter relation 
follows as a direct consequence of the consideration of the Gibbs canonical ensemble: 



which constitutes a starting point for many applications of equilibrium statistical 
mechanics. However, such a relation is only compatible with macrostates having non- 
negative heat capacities, and hence, all those macrostates with negative heat capacities 
C < cannot be appropriately described by using this statistical ensemble. In fact, such 
macrostates are thermodynamically unstable under this kind of equilibrium situations 
(a system submerged in a certain environment (heat reservoir or bath) with constant 
inverse temperature f3). 

One can easily verify from Eq.([TD that a macrostate with a negative heat capacity 
C < is thermodynamically stable provided that the correlation function {5(3ui5U) 




(2) 




(3) 
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considering the existence of correlative effects between the system and its surroundings 
obeyed the following inequality: 

(SpJU) > ks. (4) 

A simple interpretation (but not the only one possible) of the above fluctuating 
constraint follows from admitting that the thermostat or the surroundings is a finite 
system with a positive heat capacity Cui- Clearly, the existing energetic interchange 
between these systems imposes the occurrence of thermal fluctuations for the thermostat 
temperature T^, 5Ti^ = —SU/C^, with SU the amount of energy released or absorbed by 
the system around its equilibrium value. Such fluctuations can be rephrased as follows: 

Sp^ = (3HU/C^, (5) 

where the condition of thermal equilibrium /3 = /5a; is considered. By substituting Eq.Q 
into the fluctuation-dissipation relation ([1]), we obtain: 

CCro 



{5U^) . (6) 



> 1 ^ < < |C| (7) 



Finally, it is possible to arrive at the following inequalities: 

C 

by combining Eqs.(jl])-([6]). Essentially, this last result is the same constraint derived by 
Thirring in order to ensure the thermodynamic stability of macrostates with a negative 
heat capacity [9]. 

2.2. Extension of canonical Monte Carlo methods 

The study of macrostates with negative heat capacities demands that such macrostates 
be found in a stable equilibrium situation. As already discussed, such an aim could 
be implemented by considering an equilibrium situation in which the system is found 
in thermal contact with a bath having a positive and finite heat capacity C^j that 
obeys Thirring's constraint ([7]). The equilibrium condition associated with the Gibbs 
canonical ensemble is unsuitable here, since the invariability of the Gibbs thermostat 
temperature presupposes a system with an infinite heat capacity, Cuj +oo, which is 
incompatible with inequality ([7]). 

The differences between these equilibrium situations are schematically illustrated 
in FIGiH Here, the thick solid line represents the typical microcanonical caloric curve 
P (e) = ds (e) /de of a finite short-range interacting system undergoing a first-order 
phase transition, which is characterized by the existence of a regime with negative heat 
capacities (the branch p-q), with e = U/N being the energy per particle. The thin solid 
lines (e) and (e) are respectively the inverse temperature dependencies on the 
system energy per particle e of a Gibbs thermostat (with C^j +oo) and a thermostat 
having a positive finite heat capacity < C^; < +oo, with pi (e) and p2 (e) being the 
corresponding energy distribution functions. 
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Figure 1. Schematic beliavior of the microcanonical caloric /3(e) — ds{e) /de oi a 
finite short-range interacting system undergoing a first-order phase transition. Here, 
pi (s) and p2 (e) respectively represent the energy distribution functions when this 
system is placed in thermal contact with a Gibbs thermostat with inverse temperature 
i^i (e) = const and a heat bath having a finite positive heat capacity and, therefore, a 
variable (fluctuating) inverse temperature (e). 

The intersection points derived from the condition of thermal equihbrium f3 {e) = 
j3i^ (e) determine the positions of the energy distribution function p (e) maxima and 
minima. Clearly, the thermal contact with a Gibbs thermostat ensures the existence 
of only one intersection point, or equivalently, a unique peak of the canonical energy 
distribution function pi (e) for most of the admissible values of the thermostat inverse 
temperature. The important exception takes place in the inverse temperature interval 
[Pp, f3q], where there are three intersection points (two maxima Ea and Ec with C > 
and one minimum Ej, with C < 0), a fact that leads to a bimodal character for the 
distribution function pi (e) . Since no single peak can be located within the branch with 
negative heat capacities C < 0, such macrostates are poorly accessed within the Gibbs 
canonical ensemble. In fact, they turn practically inaccessible when the system size is 
sufficiently large. The existence of such a hidden energetic region constitutes the origin 
of the latent heat qi necessary for the conversion of one phase into the other during the 
coexistence of low and high energy phases (lep and hep), which are represented here by 
the coexisting peaks of the canonical energy distribution function pi [e) . 

The replacement of the Gibbs thermostat by a thermostat having a finite positive 
heat capacity crucially modifies the fiuctuating behavior and the thermal stability 
conditions of the system. In fact, one can ensure the existence of only one intersection 
point, regardless the positive or negative character of its heat capacity C, by choosing 
the appropriate thermostat and its internal conditions. In particular, it is necessary to 
ensure the applicability of Thirring's constraint ([7]) for macrostates with negative heat 
capacities C < 0. 
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The above ideas have a significant impact in the framework of Monte Carlo 
simulations. As has been discussed elsewhere |T^, large-scale Monte Carlo simulations 
are often plagued by slow sampling problems, which manifest themselves as a rapid 
increase in the dynamic relaxation time r with the system size N, causing large-size 
simulations to converge extremely slowly. These sampling problems are especially 
severe in systems near the critical point, where it is possible to distinguish two kinds 
of dynamical anomalies: (1) the so-called critical slowing down, where the relaxation 
time shows a power-law dependence on A^, r oc N"', which can be associated with the 
occurrence of a continuous (second-order) phase transition; and (2) the so-called super- 
critical slowing down, where the dynamic relaxation time exhibits a worse divergence 
with the system size: an exponential increasing r oc exp(aA^), whose incidence is 
associated with discontinuous (first-order) phase transitions. 

In Monte Carlo simulations based on a consideration of the Gibbs canonical 
ensemble, the origin of the super-critical slowing down is closely related to the existence 
of a multimodal character of the energy distribution function. Indeed, this phenomenon 
manifests itself as an effective trapping of the system macrostates in one of the coexisting 
peaks of the energy distribution function. As the system size increases, the mathematical 
form of these peaks is almost a Gaussian distribution: 



occurrence of a large energy fluctuation, whose probability p exponentially decreases as 
the system size decreases: p oc exp(— aA^). Consequently, the characteristic timescale 
for the occurrence of such rare events grows as r oc 1/p ~ exp (aA^), which explains 
the slow relaxation observed for canonical expectation values in large-scale Monte Carlo 
simulations. The existence of the above slow relaxation can be avoided if one could 
eliminate the multimodal character of the energy distribution function by appealing to 
a better control of the energy fluctuations. Fortunately, such an aim is easily achieved 
by considering a thermostat having a finite positive heat capacity C^^. 

Under this later equilibrium situation, the thermostat inverse temperature and 
the system energy (per particle) e undergo thermal fluctuations around their equilibrium 
values {(3lu) and (e), which provide a suitable estimation of the intersection point of 
the system microcanonical caloric curve (3 (e) = ds [e) / de derived from the thermal 
equilibrium condition {3 {e) = 13^ (e). Moreover, the study of the fluctuating behavior in 
terms of correlation functions {Sf3i^Se) and {Se"^) allows us to obtain the heat capacity 
C via the fluctuation-dissipation relation ([1]). Once the microcanonical caloric curve 
(3 (e) = ds (e) / de has been obtained, one can easily derive other thermodynamic 
potentials by using known integration formulae, e.g.: the entropy s {e): 




(8) 



whose width behaves 




The transition to any other peak demands the 




(9) 
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the Helmholtz free energy / (/?) = —TlogZ {f3) /N: 

Z{p)=N j exp {-N [(3e - s {e)]} de, (10) 

and the canonical averages of a certain observable O (e) : 

(O) = J O (e) exp {-N [Pe - s (e)]} de. (11) 

The simplest and most general way to implement the use of a thermostat having 
a finite heat capacity in a classical Monte Carlo calculation is through the known 
Metropolis importance sample [30]. Its extension is achieved by replacing the use of 
a constant inverse temperature in the acceptance probability: 

p{U\U + AU) =mm {exp {-(3b ^U) ,1} (12) 

with a variable inverse temperature, Pb Pu i^)- This kind of procedure can also 
be used to extend some other classical Monte Carlo methods, such as the known 
Swendsen-Wang (SW) clusters algorithm [3ll [32], [33], applicable to the Ising model 
and its generalization, the g-state Potts model: 

H,= J2 (13) 

ijSin.n 

(where n.n represents a set of nearest-neighbor lattice sites, o"j = [1, 2, . . . , g]), which 
exhibits a regime with negative heat capacities when the number of spin states q is 
greater than a certain critical value depending on the lattice dimensionality D, e.g. 
g > 3 with D=2. A direct demonstration of the applicability of the extended SW 
method using the present ideas in order to study the anomalous regime with C < in 
the 2D 10-state Potts model is shown in FIG|2l whose decorrelation time r shows a weak 
power-law dependence r ~ A^" with a ~ 0.2 at the critical point of the discontinuous 
phase transition f3c- 

Generally speaking, the consideration of a finite thermostat in order to capture the 
anomalous regime with negative heat capacities and to avoid the super-critical slowing 
down should not depend on the classical or quantum nature of the system under analysis. 
Consequently, one can expect that this idea could be used for enhancing the potentialities 
of some known quantum Monte Carlo methods. 

2.3. Complementarity character between energy and temperature 

The fluctuation-dissipation relation 1^ constitutes a particular case of a very general 
fluctuation relation: 

{Sr]SU) = kB (14) 

involving the inverse temperature difference between the surroundings (heat reservoir 
or bath) and the system rj = — (3. In fact, Eq.([l]) is obtained after substituting the 
flrst-order approximation: 

6p ~ -P^6U/C (15) 
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Figure 2. Microcanonical caloric curves of the 2D 10-states Potts raodel on a square 
lattice Lx L obtained from Monte Carlo simulations using the extended version of the 
Swendsen-Wang clusters algorithm (extended SW) . Inset panel: Power-law dependence 
of the decorrelation time r with the system size iV = at the critical point (3c, t N", 
with a ~ 0.2. 



into Eq.(fTiD. 

Alternatively, one can consider the known Schwartz inequality: 

{6A6Bf < {6A^) {6B^) (16) 

in order to rewrite the fluctuation relation (fTi|) as follows: 

A7]AU > kB, (17) 

where Ax = a/ {6x'^) denotes the thermal uncertainty of a physical observable x. Clearly, 
Eq. ffTTl) is a thermo-statistic analogy of the quantum mechanics uncertainty relation: 

AqAp > h (18) 

between position q and momentum p, which suggests the existence of certain 
complementary character between thermodynamic quantities of energy and (inverse) 
temperature [191 EQl EH ESI El] • 

It is well-known that the nature of the temperature is radically different from a 
direct observable quantity such as energy. In fact, it is a thermodynamic quantity 
whose physical meaning can only be attributed by appealing to the concept of statistical 
ensemble. In practice, the system temperature is indirectly measured by using the 
temperature of a second system through the thermal equilibrium condition, which 
plays the role of a measuring apparatus (thermometer), whose internal temperature 
dependence on some direct thermometric quantity (e.g.: electric signal, force, volume, 
etc.) is previously known. As expected, such a measuring process unavoidably involves 
a perturbation on the internal state of the system under analysis. 
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According to uncertainty relation (fTTI) . it is impossible to simultaneously reduce 
the thermal uncertainties of the inverse temperature difference Arj and the system 
energy AU to zero: any attempt to reduce the perturbation of the system energy to 
zero, AU — ^ 0, leads to a divergence of the inverse temperature difference uncertainty 
Arj — > oo, and vice- versa. Consequently, it is impossible to simultaneously determine of 
the energy and inverse temperature of a given system using the standard experimental 
procedures based on the thermal equilibrium with a second system. 

Clearly, we have to admit non-vanishing thermal uncertainties Arj and AU during 
any practical determination of the energy-temperature dependence of a given system, 
that is, its caloric curve. While such thermal uncertainties are unimportant during the 
study of large thermodynamic systems, they actually impose a fundamental limitation to 
the practical utility of thermodynamic concepts such as temperature and heat capacity 
in systems with few constituents. In order to avoid any misunderstanding, it must be 
clarified that one can obtain the energy dependence of the inverse temperature of a 
given system by calculating its Boltzmann's entropy: 

S = kehgW ^ (3 = 1/T = dS/dU, (19) 

which is possible to be achieved regardless of the system size. The limitation associated 
with the uncertainty relation ( fTTI) refers to the precision of an experimental measuring 
of the microcanonical caloric curve of a thermodynamic system. 



3. Geometrical aspects in fluctuation theory 

3.1. Starting considerations 

As previously discussed in detail in our first paper on this subject {Tj, the rigorous 
fluctuation relation (fT4l) is derived from the following ansatz for the energy distribution 
function: 

dp = p (U) dU = uj (U) n (U) dU, (20) 

where Q {U) is the state density of the system and u {U) is the probabilistic weight 
considering the thermodynamic influence of the surroundings (thermostat). Such 
functions are defined on a certain subset of Euclidean real space i?, C R : 

U G [f/jn/, Usup\- 

The next important consideration is the definition of the effective inverse 
temperature of the surroundings as: 

&(C/) = -*B^^^. (21) 

The latter assumption is not arbitrary, since it reduces to the conventional interpretation 
of this concept when one considers a closed system composed of two separable short- 
range interacting systems in thermal contact and a final thermodynamic equilibrium, 
which allows us to express the probabilistic weight uj {U) in terms of the state density 
of the second system Qb {Ub), uj ifJ) qcQb {Ut — U). 
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However, the probabilistic weight uj (U) in Eg. (1201) also admits more general 
system-surroundings equilibrium situations considering other modifying conditions, such 
as the existence of nonlinear effects driving the system-surroundings thermodynamic 
interaction, e.g., the presence of long-range interactions [3ll[35], or a system acting as 
the surroundings that is found in a metastable equilibrium whose relaxation time is 
so long that its dynamic evolution can practically be disregarded, such as the case of 
systems with glassy dynamics [36j. Such circumstances explain why we refer to the 
inverse temperature ( l2Ti) as effective. 

The number of microstates W used to obtain the Boltzmann's entropy (flQl) is given 
by the coarsed grained definition: 

W = nSe, (22) 

with 6e being a certain small constant energy that makes W dimensionless. The work 
hypothesis ( l20l) and definitions (12T|) and (l22l) allow us to express the inverse temperature 
difference as: 

V (U) = (3^ (U) ~(3{U)^ -kB^-^^^l^- (23) 
With the above relation, one can easily obtain the thermodynamic identities: 

fU sup 

{V)= / viU)piU)dU = 0, (24) 

r Usup 

{Ut]) = / Ur,{U)p{U)dU = kB, (25) 

which are derived by integrating by parts and considering the following boundary 
conditions: 



p (Uinf) = p (Usup) = (26) 

0. (27) 



dp{Uinf) _ dp (Usup) 



dU dU 

Eq. (1241) is simply the thermal equilibrium condition expressed in terms of statistical 
expectation values: 

w = m . (28) 

This rigorous result clarifies that the known equalization of (inverse) temperatures 
during the thermodynamic equilibrium of two systems actually takes place in an average 
sense. The fluctuation relation of Eq. flT^ is obtained from Eqs. fl2^ and fl2^ after using 
the identity {5U5r]) = {Ut]) - {U) {r]). 

Let A be a continuous and differentiable function on M^, which also admits a 
bound expectation value (A), \{A)\ < +oo. Under these assumptions, one can obtain 
the following thermodynamic identity: 

ks—\ = {Ar^)^{5A5r^). (29) 
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In particular, this identity reduces to Eg. (121]) and Eg. (1251) for A = 1 and A = U 
respectively. Moreover, it also drops to the remarkable fluctuation relation: 

'f^B^ + ^V') = (30) 



'dU 

for A = rj. This latter identity, hereafter referred to as the complementary fluctuation 
relation, accounts for the same information about the system stability conditions 
derived from the fluctuation relation of Eg. (fT^ . For instance, by using the Gaussian 
approximation (see subsection 13.41 below): 

■dr,{U)\^dr,[U) 

'dir/-~dir ^^^^ 

and focusing on the eguilibrium situation between two separable short-range interacting 
systems, we obtain the fluctuation relation: 

^"^^^B = , (32) 
which leads to the same stability criterion derived from Eg. ([6]). 



3.2. Reparametrization invariance 

Let us consider a bijective application G (U) : R —>■ R, which is a piece-wise continuous 
and two time differentiable function of the variable U. Such a function allows for 
the existence of a bijective map : — > of the subset on another subset 

C -R : 6 G [Qinf, Qsup]- It could be said that these subsets constitute two eguivalent 
coordinate representations of all admissible macrostates of the system, which shall be 
denoted as Ru and R^ respectively. The coordinate transformation induced by the 
bijective function B (U) is referred to as a reparametrization. 

Since the elementary subset [U, U + dU] represents the same system macrostates 
considered by the elementary subset + dQ], the elementary probability dp that 

the system is found in such conditions, Eg. (l2Ul) . does not depend on the coordinate 
representation used for its expression: 

dp = p„ (U) dU = (6) de. (33) 

Here, pu (U) and (0) denote the system distribution functions in the representation 
Ru and R^, respectively, which are mutually related by the transformation rule: 

(34) 

Let dW be the elementary volume considering the number of microstates belonging 
to the elementary subset [U, U + dU]. As the case of the elementary probability dp, dW 
does not depend on the coordinate representation, and hence, it obeys the following 
properties: 



P0 (6) = p„ (U) 



dQ 
dU 



dW = Qu (U) dU = (6) dQ, 



(35) 
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(6) = {u) 



dU 



(36) 



Consequently, the probabilistic weight Uu iU) considering the surroundings thermody- 
namic influence behaves as a scalar function under reparametrizations: 

Uu{U)=uj^{Q)=uj^[U{Q)], (37) 

with U (6) being the inverse function of 6 {U). 

The reparametrization invariance of the probabihty distribution function also leads 
to the reparametrization invariance of the expectation value of any physical observable 
A = A{U) = A{Q) (scalar function): 



A{Q)p^{e)de 



in! 



A{U)p^{U)dU={A)^. 



{31 



(39) 



Uinf 



such that, one can denote the expectation values without indicating the coordinate 
representation used for its expression: 

= {A)^ ^ {A) . (40) 

A remarkable equilibrium situation of the conventional thermodynamics and 
statistical mechanics is the system in energetic isolation, whose probabilistic weight: 

1 



{U\Uq 



-6{U-Uo) 



(41) 



defines the known microcanonical ensemble. This probabilistic weight possesses the 
notable feature that its mathematical form does not depend on the representation: 

1 



{U\Uo 



(eieo) 



-5(0-60 



^^<^ (Oo) 

a property that is straightforwardly derived from the identity: 

56 (t/o) 



5(0-eo) =(5(t/-f/o) 
and the transformation rule (l36l). 



(42) 



(43) 



3.3. Reparametrization duality 

Let us define the thermostat inverse temperature in representation as: 

/3^ 



dloguj^, (6) 



de 

Therefore, it obeys the transformation rule: 



dU 



(44) 



(45) 



as a consequence of the scalar character of the probabilistic weight u^, with being 
the thermostat (effective) inverse temperature expressed in Eq. fl21l) . 
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Boltzmann's entropy of the system in this representation can be defined by: 

S^ = kBlogW4„ (46) 

where Wip = Q^p^e^, with ^e^ being a suitable constant that makes Wip dimensionless. 
The above coarsed-grained definition of Boltzmann's entropy is not properly a scalar 
function as the case of the probabilistic weight ( !37|) . In fact, it obeys the transformation 
rule: 

f dQ 5eu \ / . ^,^ 

' ' (47) 

the density 
(48) 



S^ = Su- ks log 

As already pointed out by Ruppeiner (see subsection II. B of ref. [25j' 
distribution function derived from Einstein's postulate: 

-S{x)- 



Px{x)dx = C exp 



kf 



dx 



obeys different mathematical forms under different coordinate representations, x — >■ 
y{x), if one assumes that the entropy is a state function whose value does not depend 
on the representation (scalar function), S{x) = S{y). A simple analysis allows us to 
verify that the left-hand side of Eq. ( l48l) behaves as: 

dx{y) 



Px{x)dx = p{x) 

while its right-hand side as: 
-S{x) 



dy 



dy = py{y)dy, 



C exp 



k 



B 



dx = C exp 



S{x) 



k 



B 



dx{y) 



dy 



dy 7^ C' exp 



S{y) 



iB 



(49) 



dy. (50) 



This fact not only constitutes an important defect in order to develop a Riemannian 
formulation of fiuctuation theory, but it also presupposes some inconsistences with the 
thermodynamic arguments behind of Einstein's postulate for the fiuctuation formula of 
Eq. fHHl) . In this work, we shall assume the entropy modification fH7|) associated with 
reparametrizations and analyze its direct consequences. Clearly, such an alternative 
definition allows us to preserve the functional dependence of fiuctuation formula (l48ll 
in any coordinate representation. It requires that the entropy is no longer a state 
function with a scalar character, as is usually assumed in other geometric formulations 
of fiuctuation theory [25] . 

Under these above assumptions, the system inverse temperature /J"^ in the 
representation is given by: 



and obeys the transformation rule: 



d_ 

'dU 



(dQ 5eu\ 




'dQ' 


\dU5eJ_ 




_dU_ 



(51) 



(52) 



As expected, the inverse temperature difference in the representation R^ can be 
expressed as: 



(0) = Pt (0) - P"" (0) = -k 



B- 



(53) 
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By only admitting regular reparametrizations obeying the constraints: 



< 



dU 



< +00, < 



< +00 



(54) 



on every point U G M„, one can easily show the validity of the boundary conditions: 



P<t> (©in/) — P<t> (©sup) — 0, 



<9p0 (© 



inf) 



sup J 







(55) 
(56) 



de de 

by starting from Eg. (1261) and Eq. (|271) . 

As already shown in the previous subsection, definition fl53|) and the boundary 
conditions (l55l) and ( !55|) lead to the following extensions of the rigorous identities (l24l) . 
((25]) and ([29]): 

{V<p)= r'"%40) p40) rf© = 0, (57) 



©in/ 



; A 



as well as the generalized fluctuation theorems: 



-k 



B 



dQ 



0, 



(58) 
(59) 

(60) 
(61) 



and finally, the thermodynamic uncertainty relation: 

A0Ar/^ > kB. (62) 

Thus, the consideration of coordinate changes makes it possible to extend the results 
already derived by using the energy representation Although the thermodynamic 
identities fl2ll25D and fl57ll58D. and the fluctuation theorems ffT4ll30D and fl60]l6TD. as well 



as the uncertainty relations (fTTl) and (!62l) are closely related, they represent different 
thermodynamic relations characterizing the same equilibrium situation. It could be said 
that all of these mutually related identities account for the existence of a special kind 
of internal symmetry, which shall be hereafter referred to as reparametrization duality. 

The invariance under reparametrizations (coordinate transformation or diffeomor- 
phisms) is the same kind of symmetry considered by Einstein's theory of gravitation. 
However, there exist radical differences between this latter physical theory and the geo- 
metric statistical formalism developed in this work. (1) While the gravitation theory is 
defined in terms of local quantities, the rigorous thermodynamic identities obtained here 
are expressed in terms of statistical expectation values defined over the entire subset 
representing all admissible system macrostates in the present equilibrium situation, that 
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is, this is a non-local theory, similar to quantum mechanics. (2) Furthermore, Einstein's 
theory refers to the same physical laws in different representations, while the above 
thermodynamic identities consider a family of different fluctuations relations exhibiting 
the same mathematical appearance under different coordinate representations of a given 
equilibrium situation. This is why we refer to it as reparametrization duality instead of 
reparametrization symmetry. 

In the next subsection, we shall arrive at a local formulation of the present approach 
with a Riemannian-like structure closely related to other geometric approaches of 
fluctuation theory existing in the literature [25]. We shall see, however, that such a 
development presupposes the consideration of certain unexpected approximations. 

3.4- Riemannian approach 

Let us assume that the systems under consideration are large enough to deal with 
the thermodynamic fluctuations by using a Gaussian approximation. An essential 
assumption considered here is that the system undergoes small thermal fluctuations 
close to its equilibrium point, which is determined by the most likely macrostate. 

A problem encountered is that the most likely macrostate actually depends on the 
coordinate representation used for describing the system behavior, which is a direct 
consequence of the non-scalar character of the system entropy. In order to show this 
fact, let us consider the transformation rule of the inverse temperature difference: 



The stationary condition associated with the most likely macrostate in each 
representation are given by: 



According to Eq. (!63l) . the vanishing of rju does not correspond to a vanishing of rj^, and 
vice versa, a result showing that the most likely macrostate depends on the coordinate 
representation. 

This last result contracts with the general validity of the thermal equilibrium 
condition in terms of statistical expectation values, Eq. flFTI) . It clearly indicates that 
the method generally used for deriving such a condition in terms of the most likely 
macrostate is just a suitable approximation. Nevertheless, it could be easily noticed 
that the modification involved during the reparametrization change is just a second- 
order effect. The transformation rule ( !63l) can be combined with Eq.( l24l) and Eq.( !57l) 
in order to obtain: 




(63) 



r]u{U) =0 for Ru and r]^ (©) = for R^, 



(64) 




(65) 



where the following notation is considered: 



_ de iu) 



dU 



(66) 
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Eq. (1651) indicates that the second additive term on right-hand side of the transformation 
rule fl63l) is just a small correction, which can be disregarded in most practical 
applications. Therefore, one can admit the approximate relation: 

r^^ = r^„(At)-'^0, (67) 

where A denotes the value of the function A (U) at the most hkely macrostate, 
A = A{U). 

Basically, the approximation assumed in Eq. fl67|) is equivalent to considering 
Boltzmann's entropy ( H6l) as a scalar function, and hence, the approximate 
transformation rule of the system inverse temperature is given by: 

= {Aty\ (68) 

In general, the Gaussian approximation allows us to consider the fluctuations of an 
arbitrary energy function A (U) as: 

6 A = —^^U- (69) 
In particular, it allows us to introduce the following transformation rule: 

56 = Kt5U. (70) 
Moreover, by starting from Eq. (!67|) . we obtain: 



= [^Vu - Vu^ log [kt^-^^ SU 



(71) 



which reduces to: 

= {Aty' 6vu, (72) 

after considering the thermal equilibrium condition f]^ = 0. Using this latter 
transformation rules, one can obtain the transformation rules of some fluctuations 
relations: 

{6e') = {AtY {6U') (73) 

{H)=C^tnSvl) (74) 
{6067]^) = {6U67]^) = kB. (75) 

Exactly, Eqs. fl73l) - fl75l) correspond to transformation rules of contravariant second- 
rank tensors, covariant second-range tensor and scalar functions in a differential 
geometric theory, respectively. In order to provide a Riemannian structure to the present 
geometrical approach, we must introduce an appropriate metric. Such a role could be 
carried out by the global curvature K^: 

dr]^ d'^ log P4, 

= ae = -^^^e^ ^^^^ 
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evaluated at the most likely macrostate, which allows for the conversion between the 
fluctuations of the conjugated thermodynamic quantities (covariant and contravariant 
vectors) within the Gaussian approximation: 

6r]^ = K^6e. (77) 

The global curvature obeys the transformation rule: 

K^ = {At) Ik., 



dc4, 



dU 



with C(f) = log [A^Seu/Setp) , which reduces to: 

K^={nY^K^ (79) 

after considering the thermal equilibrium condition f]^ = ^ and dismissing small 
contributions associated with the non-scalar character of Boltzmann's entropy (the two 
terms associated with the Boltzmann's constant fc^). Clearly, the global curvature can 
only be considered as a second-rank covariant tensor under the above approximations, 
since the general transformation rule (I78|) does not correspond to this kind of geometric 
object. Interestingly, such a function appears in the complementary fluctuation relation 
(!6T|) . which establishes the non-negative character of its expectation value in any 
coordinate representation: 

kB {K^) = {6vl) . (80) 

As already commented, this rigorous fluctuation relation satisfies, as a whole, the 
reparametrization duality, which is not the case of the global curvature considered 
as an individual entity. 

By using the global curvature K^, one can easily obtain other fluctuations relations 
such as: 

{Sr]t5Q) = (502) = {5U^) = {5r^l5U) = kj,, (81) 
and rewrite the distribution function in this Gaussian approximation as follows: 



(e|e) dQ 



27rk, 



■ exp 



2k, 



de. 



12) 



4. Generalized Gibbs canonical ensemble 



Let us denote by the thermostat temperature in the representation R^, with 
= l/T^. One can formally introduce the heat capacity of this representation 

as: 



(83) 



which allows us to obtain a geometric extension of fluctuation-dissipation relation ([T]): 



kBC^ = {Mf {5Q^) + C^{5(5t5Q) 
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after combining the Gaussian approximation: 

= -Pl/C^SQ (85) 

with definition ( l53l) and the fiuctuation relation ( l60l) . A relevant case among 
the admissible equilibrium situations considered by the above fiuctuation-dissipation 
relation is the one obeying the constraint 6f3^ = 0, which is associated with the following 
distribution function: 

dPc ) = -^-X (-^f^^^) (®) (86) 

This is just the analogous version of the Gibbs canonical ensemble in the 
representation, with f]f being a constant parameter. By rewriting this particular 
distribution function in the energy representation i?„: 



dPc {U\Pt)=^-^exp 

z ( /3r 



-^(3tQ {U) 



(U) dU (87) 



one arrives at the same expression found for the so-called generalized canonical ensemble 
recently proposed in the literature [261 EZ] • Let us now analyze its general mathematical 
properties. 

4.I. General mathematical properties 

As usual, the partition function Z derived from the normalization condition: 

Z (pf) = r^^' e"^^'®^^^fi„ (U) dU (88) 

allows us to obtain the generalized Planck's thermodynamic potential: 

P^{Pt)=-kBlogZ{Pt), (89) 
which provides two relevant statistical expectation values: 

These last results can be combined in order to obtain the canonical version of the 
fiuctuation-dissipation relation (1841) : 



- ^^^-3 = (^0'> ^ ^B^l = ^Ptf {se') , (91) 

with being the canonical heat capacity: 

- dT^ ■ ^^^^ 

Clearly, this theorem states that the stable thermodynamically macrostates are those 
with a nonnegative heat capacity > 0. 
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Let us now rewrite Planck's thermodynamic potential in the representation: 

g-P,(/3c^)Afl ^ /"°'"'e~^^'®^]^(e)c/0 (93) 



and develop a Gaussian approximation (the second-order power expansion in 0) around 
the local maxima: 



P 



dQ 



e-m^- I ' ^-^^'^^"""^ (95) 

with AG = — 6c and given by: 

P; = mi{PtQ-S^{Q)]. (96) 

The local maxima 6c are derived from the stationary and stability conditions: 

f^c = — QQ — = i^c) , = — > 0. (97) 

By admitting the existence of only one maximum, this approximation yields: 

1 , / 2Tiki 



p,m-^p;-l^og^], (98) 



(A62> = ks^. (99) 

Clearly, the additive logarithmic term in the Gaussian estimation of the Planck 
thermodynamic potential constitutes a small correction in the case of sufficiently large 
systems. By dismissing this small contribution, one finds that Planck's thermodynamic 
potential is approximately given by the known Legendre transformation: 

p;{Pt) = mf{pfe-s,{e)}. (loo) 

The stationary condition is merely the condition of thermal equilibrium associated with 
this representation, while the stability condition is simply the requirement of non- 
negativity of the microcanonical heat capacity C^: 

^^^ = -(^*yj;<o^c,>o. (101) 

Eqs.( l88l) -( 110ip correspond to many well-known dual expressions previously obtained 
within the Gibbs canonical ensemble ([3]) . Obviously, these two ensembles are intimately 
related. By considering the scalar character of the probabilistic weight ujf 

(6|/?c") = -^-r exp (-^PtQ) , (102) 

the thermostat inverse temperature (3^ in the energy representation Ru is given by: 

(^:iU)=pf-^^- (103) 
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This latter result clarifies that the generalized canonical ensemble fl87p corresponds to 
a special kind of equilibrium situation with a variable (fluctuating) inverse temperature 
of all admissible states accounted for by fluctuation-dissipation relation ([T]), that is, a 
situation with non-vanishing system-surroundings correlative effects {6f3^6U) ^ 0. 
By considering the transformation rule for the microcanonical curvature k^: 

d^c^ _ / ^ 
dU ' '^^ dU^ [dU 



(104) 



one can flnd that the requirement > can be combined with the existence of 
macrostates with k„ < in the energy representation with an appropriate selection 
of the reparametrization U ^ Q {uj^ . This fact is more evident when working in the 
energy representation Ru, where the stability condition reads as follows: 

+ i^l^^ = -^u^ > 0. (105) 

By considering the relations = {^""fjCu and d^HdV = (/3^)VC with C„ and Q 
being the heat capacities of the system and the thermostat respectively (their usual 
deflnitions), as well as by using the thermal equilibrium condition = = /3, one 
arrives at the expression: 

>0, (106) 



which leads to Thirring's stability condition ([7]) for macrostates with Cu < 0. 

As the Gibbs canonical ensemble (I3l), the present geometric extension ( IHTl) becomes 
equivalent to the microcanonical ensemble with increasing of the system size N, 
A0/O ~ an equivalency that can be ensured even for macrostates with 

Cu < or K„ < with an appropriate selection of the reparametrization 0(f/). 
This remarkable property makes this ensemble a very attractive thermo-statistical 
framework, since besides of exhibiting many notable properties of the usual the Gibbs 
canonical ensemble, it also provides a better treatment of the phenomenon of ensemble 
inequivalence associated with the presence of negative heat capacities, as already 
discussed in refs. [261 122] • In particular, this statistical ensemble constitutes a suitable 
framework for extending of Monte Carlo methods, as discussed in subsection I2.2[ 



4-2. Derivation from information theory 

It is possible to realize that the generalized Gibbs canonical ensemble ( ISTl) can also 
be derived from Jaynes's reinterpretation of statistical mechanics in terms of the 
information theory of Shannon [28], e.g., by considering the maximization of the known 
statistical (extensive) information entropy: 

Se = - ^Pk^ogpk, (107) 

k 

I The presence of additive terms with Bohzmann's factor /cs in Eq. (|104p takes into account the 
modification of the system entropy during a reparametrization and the consequent correction of the 
most hkely macrostate. 
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under the normalization condition: 



(1) = Ep'^ = 1 



(108) 



k 

and the following nonlinear energy-like constraint: 



(e) = ^e(f/,)p,. 



(109) 



k 



Such a derivation was developed by Toral in ref . [27] . The interested reader can refer to 
this work for more details. 

Clearly, the bijective character of the reparametrization U ^ <d (U) should ensure 
that this generalized ensemble exhibits almost the same stationary properties obtained 
from the application of the Gibbs canonical ensemble in sufficiently large systems, where 
one usually assumes the appropriateness of the Gaussian approximation. However, the 
nonlinear character of the bijective application (U) produces a deformation in the 
canonical description, which conveniently modifies the system fiuctuating behavior and 
the accessible regions of the subset of all admissible system macrostates M^. 

4-3. Connections with inference theory: generalization of Mandelbrot's approach 

Generally speaking, statistical inference can be described as the problem of deciding 
how well a set of outcomes {xi,X2, ■ ■ ■ ,Xm,), obtained from independent measurements, 
fits to a proposed probability distribution: 



If the probability distribution is characterized by one or more parameters {6), this 
problem is equivalent to inferring the value of the parameter (s) from the observed 
measurement outcomes x. To make inferences about the parameter, one constructs 
estimators, i.e., functions: 



of the outcomes of m independent repeated measurements [37]. The value of this 
function represents the best guess for 6. 

Commonly, there exist several criteria imposed on estimators in order to ensure 
that their values constitute good estimates of the parameter 6, such as: 

• Unhiasedness: 




(110) 



e{Xi,X2,...,Xm) 



(111) 




m 



(112) 



k=l 



• Efficiency or minimal statistical dispersion: 



/2 
{o-{d)) \{dp{x,\e) 



minimum. 



(113) 



k=l 
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• Sufficiency: 

dp {xi,X2, ...,x^\9) = f{xi,X2,..., Xm) dp (9^ , (114) 

where dp (^9^ is the marginal distribution of 9 and / {xi,X2, ■ ■ ■ , Xm) is an arbitrary 
function of the measurements, independent on 9. 

Since any statistical estimator 9 represents a stochastic quantity, it is natural in 
inference problems that an estimator obeys the unbiasedness (11121) and efficiency (11131) 
conditions. However, there exists a remarkable theorem of inference theory, the Cramer- 
Rao 's inequality, which places an inferior bound on the efficiency of an arbitrary unbiased 
estimator: 

^ i^y 

where Ip {9) is the so-called Fisher's information entropy: 



h{9) 



d\ogp {x\9) 



2 



p{x\9)dx. (116) 



89 

On the other hand, efficiency condition (I114p ensures that, given the value of 
^ (xi, X2, . . . , Xm)-, the values of the data (a:i, X2, . . . , Xm) are distributed independently of 
9, containing in this way all of the information about parameter 9 that can be obtained 
from the data. As with unbiasedness and efficiency, sufficiency is also a natural desirable 
condition in inference problems. However, a theorem by Pitman and Koopman [38] 
states that sufficient estimators only exist for a reduced family of distribution functions, 
the so-called exponential family: 

dp {x\9) = exp [A (9) + B{x)C (9) + D (x)] dx. (117) 

Mandelbrot was the first investigator to realize the intimate connection between 
statistical mechanics and inference theory [21]. Clearly, the Gibbs canonical ensemble ([3]) 
constitutes a relevant physical example of probabilistic distribution function belonging 
to the exponential family (11171) . As the well-known Kinchin work in the framework 
of information theory [39], Mandelbrot proposed a set of axioms in order to justify 
a direct derivation of the Gibbs canonical ensemble in the framework of inference 
theory. Moreover, he also focussed the inference problem of the inverse temperature 
/?, which appears as a parameter of the Gibbs canonical ensemble ([3j), through some an 
unbiased estimator /3 defined for a set of outcomes of the system energy U. Thus, this 
author provided an interpretation of the energy-temperature complementarity previously 
postulated by Bohr and Heisenberg [IHl [20] : 

AJA.U > kn, (118) 

with AcX = a/ {Sx"^)^, a result that follows from the Cramer-Rao's inequality (11151) 
after noting that the Fisher's information entropy (11161) for the Gibbs canonical 
distribution ([3]) is simply the canonical expectation value (*)^ of the energy dispersion. 
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After reading the present discussion, one can point out some critiques to 
Mandelbrot's approach. In regard to his interpretation of energy-temperature 
complementarity, Eg. (11 181) . it is clear that such an uncertainty relation only applies 
in the framework of the Gibbs canonical ensemble ([3]). Moreover, this inequality 
accounts for the limits of precision of a statistical estimation of the inverse temperature 
P appearing as a parameter of the canonical ensemble ([3]). Clearly, this quantity has 
nothing to do with the system inverse temperature, but rather the inverse temperature 
of the Gibbs thermostat. This is a common misunderstanding of some contemporary 
developments of statistical physics, where it is not distinguished between these two 
temperatures, leading in this way to some limitations and inconsistences. Clearly, such 
difficulties are overcome by the uncertainty relation f|T7|) associated with the energy- 
temperature fluctuation-dissipation relation ([T]). 

The differences between the Gibbs temperature of the canonical ensemble ([3]) 
and the Boltzmann's deflnition ([5]) are irrelevant in the case of large short-range 
thermodynamic systems considered in conventional applications of statistical mechanics 
and thermodynamics, overall, in those physical situations where the necessary conditions 
for the equivalence between canonical and microcanonical descriptions apply. However, 
the existing differences become critical when one considers the thermodynamical 
description of long-range interacting systems such as the astrophysical ones, where 
the presence of macrostates with negative heat capacities constitutes an important 
thermodynamic feature that rules their macroscopic behavior and dynamical evolution 
[T3l [H]. As already discussed, such an anomaly cannot be described by using the 
Gibbs canonical description ([3]). Besides, there does not exist in this context an 
appropriate Gibbs thermostat that ensures the existence of a thermal contact (a 
boundary interaction) in presence of a long-range interacting force such as gravity. 

The above limitations also extend to other physical contexts such as small or 
mesoscopic nuclear, molecular and atomic clusters, where the presence of a negative 
heat capacity is not an unusual feature [13], while the thermodynamic influence of a 
Gibbs thermostat constitutes a very strong perturbation of its internal thermodynamic 
state. In this kind of scenario, there does not always exist a clear justiflcation for 
the direct application of some theoretical developments based on the consideration of 
the Gibbs canonical ensemble, e.g.: the use of flnite-temperature calculations for the 
study of collisions in high energy physics. Interestingly, a collective phenomenon such 
as the nuclear multi-fragmentation resulting from collisions of heavy nuclei is simply a 
flrst-order phase transition revealing the experimental observation of macrostates with 
negative heat capacities C < P [5] . Clearly, such a realistic phenomenon cannot be 
appropriately described by using the canonical ensemble. 

Remarkably, its is easy to note that the Gibbs canonical ensemble ([3]) is not the 
only one probabilistic distribution function justifled in terms of inference theory, as 
originally presupposed by Mandelbrot in his approach. In fact, the whole family of 
the generalized Gibbs canonical ensembles ( 187|) also belongs to the exponential family 
( 11171) . and hence, such distributions also ensure the existence of sufficient estimators 
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obeying uncertainty relations a la Mandelbrot: 



A^pfA^e > kB, 



(119) 



as a consequence of the underlying reparametrization duality discussed in this work. As 



expected, A(^x = \J {^x"^) ^, with (*)^ being the generalized canonical expectation values 
derived from the generalized ensemble (1871) . 

5. Conclusions 

We have provided in this work a panoramic overview of direct implications and 
connections of the energy-temperature fluctuation-dissipation relation ([1]) with different 
challenging questions of statistical mechanics. 

As briefly discussed, the main motivation and most direct consequence of this 
generalized fluctuation relation was the compatibility with macrostates having negative 
heat capacities in the framework of fluctuation theory. Such a feature makes possible 
to analyze and apply the necessary conditions for the thermodynamical stability of 
such anomalous macrostates in order to extend the available Monte Carlo methods 
based on the consideration of the Gibbs canonical ensemble ([3]), a procedure that 
also allows one to avoid the incidence of the so-called super- critical slowing down 
encountered in large-scale simulations. Moreover, the fluctuation-dissipation relation 
constitutes a particular expression of a fluctuation relation leading to the existence of a 
complementary relationship between thermodynamic quantities of energy and (inverse) 
temperature (|T7|) . 

The consideration of geometric concepts, such as coordinate changes or 
reparametrizations, leads to a direct extension of many old and new rigorous results 
of statistical mechanics in terms of a special kind of internal symmetry that we refer to 
here as a reparametrization duality. Such a basis inspires the introduction of a geometric 
generalized version of the Gibbs canonical ensemble fl87|) . which has been recently 
proposed in the hterature [26l [27]. This latter probabilistic distribution allows for a 
better treatment of the phenomenon of ensemble inequivalence or for the consideration 
of anomalous macrostates with negative heat capacities. At the same time, this family 
of distribution functions still preserves many notable properties of the Gibbs canonical 
ensemble, including its derivation from Jaynes' reinterpretation of statistical mechanics 
in terms of information theory, as well as Mandelbrot's approach based on inference 
theory. 
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